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In the context of theories where particles can have different limiting velocities, we review the 
running of particle speeds towards a common limiting velocity at low energy. Motivated by the 
recent OPERA experimental results, we describe a model where the neutrinos would deviate from 
the common velocity by more than do other particles in the theory, because their running is slower 
due to weaker interactions. 



1. INTRODUCTION 

Recently, we [lj considered theories where the parti- 
cles have different limiting velocities. Due to interac- 
tions, these velocities become running parameters and 
we showed that the limiting velocities are driven by the 
renormalization group (RG) equations towards a univer- 
sal limiting velocity - the "speed of light" - at low energy. 
Because the running is produced by the self-interaction 
diagrams involving the other particles in the theory, we 
noted that the differences between the limiting speeds is 
greatest when the interactions are the weakest, and we 
pointed to gravitational waves as potentially the most 
sensitive test of this idea. 

Recently the OPERA collaboration has presented evi- 
dence that neutrinos may have a different limiting veloc- 
ity than the speed of light Q ■ Neutrinos are the second 
weakest interacting particles, after gravitons. While we 
are aware of potential problems with the experiment Q 
and potential difficulties for the idea of a neutrino speed 
that is different from c we would like to address 

here a mechanism for allowing the neutrinos to have a 
greater velocity difference from c than the constraints on 
other particles would otherwise appear to allow. 

The ideas of renormalization group evolution of the 
limiting speeds and the general approach towards a com- 
mon speed at low energy are independent of the under- 
lying mechanism, as long as one has an effective four- 
dimensional theory of particles with different speeds in- 
teracting with each other. However, our original moti- 
vation involved emergent fields. The Weinberg- Witten 
theorem is usually interpreted as telling us that non- 
Abelian gauge bosons and gravitons cannot be emergent 
fields arising from any underlying Lorentz invariant non- 
gauge theory. Such a Lorentz- violating underlying theory 
would generally involve wave equations which have dif- 
ferent limiting velocities, so it would be a puzzle as to 
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why Lorentz invariance appears as such a good symme- 
try at low energy. Our work using the renormalization 
group [l[ is a partial answer to the issue of an emergent 
Lorentz symmetry. 

In section 2 we briefly review the running of the lim- 
iting speeds. In section 3 we present a model in which 
the neutrino limiting speed differs from c by a greater 
amount than does that of the charged leptons. A brief 
discussion completes this note. 



2. LIMITING VELOCITIES AS RUNNING 
PARAMETERS 

In a theory with Lorentz invariance, there is a common 
limiting speed c and this speed can be used to define a 
natural set of units in which c = 1 . However if the theory 
does not have a full Lorentz invariance, the wave equa- 
tions for the various fields can involve limiting velocities 
Cj that are different from each other. Because these veloc- 
ities are no longer universal, the limiting speeds cannot 
be universally removed from the theory and differences 
among them become important. Indeed, quantum cor- 
rections may modify the propagation of all the particles 
through the self-energy diagrams. 

Self-energy diagrams generally involve more than one 
type of field. For example, in QED the photon self-energy 
involves fermions and the fermion self-energy involves a 
photon in addition to a fermion. Due to the lack of uni- 
versal Lorentz invariance in such a theory, the time-like 
component of the propagator and the space-like compo- 
nent get renormalized differently Because the limiting 
speed is related to the coefficients of energy and momen- 
tum variables, this means that the limiting speed itself 
gets renormalized, and the renormalization depends on 
the limiting velocity of both the fermion (c/) and the 
gauge boson (c s ). Lorentz invariance would guarantee 
that there is no running of these parameters in the limit 
that Cf — c g , so that we know that the running will be 
proportional to c/ - c 9 in the limit of small differences. 
The parameters do run when c/ — c g ^ 0. To decide 
whether the speeds run towards each other, or away from 
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each other, at low energy require s explicit calculations. 

All of the cases that we treated in [l| resulted in the 
speeds approaching a common velocity at low energy. For 
example, the case of an interacting fermion and U(l) 
gauge boson resulted in the beta functions 
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c g as an infrared attractive fixed line. 



These lead to c/ 

The rate of approach towards the fixed point depends 
on the strength of the interaction and the number of fields 
involved. For a single field, the approach is logarithmic, 
which raises the issue of whether approach is sufficient 
to satisfy experimental constraints. In [lj we provided a 
model with more fields, elaborated below, in which the 
running was much faster. The issue of whether the run- 
ning would be logarithmic over large energy scales, or 
faster over a narrower energy range, is model dependent. 



3. THE MODEL 

One of the proposals to account for superluminal neu- 
trinos is to use Coleman- Glashow idea [8j that any par- 
ticle will have its own mass as well as maximum at- 
tainable velocity. Although this scenario is in perfect 
agreement with the long baseline neutrino experiments 
with Vu/cry « 1 + 10~ 5 , where the best fit with the data 
leads to energy independent profile [j| , it does not tell us 
why there is a stringent bound on the difference between 
the maximum attainable velocity of charged leptons and 
speed of light @ 
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Here we argue that if Lorentz symmetry is an emergent 
symmetry with a large number of fields, one can explain 
the hierarchy between the limiting speeds of electrons 
and neutrinos. 

One of the keys to this puzzle is to notice that neu- 
trinos, unlike all other leptons, are not charged under 
^em(l)- Now, let us introduce a large number N 7 of hid- 
den U(l) photons as well as a large number Nz copies of 
neutral Zs (under the U(l)s) in addition to the Standard 
Model (SM) ones [10J. Moreover, we assume that there 
are a few leptons Ni v (hidden and SM) that we assume 
for simplicity have the same origin, and hence all have 
the same initial speed of light 1 + A, with \cj\ << 1, at 
some UV emergent scale fx*. The charged leptons, de- 
noted by £, have a common initial charge Z* a under the 
different E/(l)'s gauge sectors, while both charged and 
neutral leptons (the neutral ones are denoted by v) are 
charged initially with charge Z\ v z under the Nz copies 
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fermions. At the UV scale the fermions and gauge bosons 
are taken to be massless and hence will participate in the 
RG running. As we run down our RG equations, most of 
the hidden particles become massive and decouple from 
the RG system. This situation is illustrated in FIG. [U 
Here we do not assume any specific model for the nature 
of the coupling between the leptons and Zs. All we need 
to assume is the existence of a positive sign /3 function 
for f3(Z^ uZ ). In order to produce the required hierarchy 
we take A 7 >> N z » N e . 







FIG. 1: Lower order quantum loops for iV 7 photons and Nz 
neutral Z's interacting with Ne leptons. 

Since all fermions share a common initial speed (also 
all gauge bosons have their common initial speed) and 
common initial coupling strengths, the evolution of the 
system can be modeled with only six parameters: c 7 , 
cz, C£, and c v are respectively the photon, Z, charged 
lepton, and neutrino speeds, in addition to the coupling 
strengths Zu n and Zg v ^z- Writing the speed of species i 
as Cj = 1 + 5i, we find that the (RG) equations read 
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where /3<s 7 = (4-n-) 2 fidS^ / dfx ate, and {gi} are 0(1) loop- 
numbers that can be calculated in any realistic setup. 

As was shown in [l|, the positivity of {<?;} suffice for 
the existence of a stable IR fixed line c 7 = cz = ce = c„. 
In the following we perform a series of approximations 
that enable us to show our main point. Namely, the 
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hierarchy 7V 7 >> Nz » N& leads to a hierarchal struc- 
ture of the different speeds with respect to the speed 
of light. Some of these approximations will turn out to 
be crude when compared to numerical integration of the 
whole system. However, this procedure is sufficient to 
illustrate our point. Subtracting Eqs. [3] and 0] and El 
and 0] and [51 and using N y » Nz >> Nt, therefore 
neglecting the less dominating terms, we obtain 
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Finally, we compute j3 (<5„ — <5 7 ) by noticing 

P(S V -S 7 ) = f3{S„-Sz) + P(S z -S e ) 
+/3 (S e -S 7 ) 
w g 6 N z Z^ z (Su - Sz) 
-g 5 N z Zj z ( s t - Sz) ■ 
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Eqs. O to [T2] as well as [7] and [5] can be integrated directly 
to obtain 
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where, for example, 
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It turns out that the corrections in the above expression 
are important, despite being subleading, within the nu- 
merical simulations. However, this shows that charged 
leptons, that see iV 7 , have more enhanced running than 
neutrinos, which only see N z << N 7 . 

In FIG.[3]we demonstrate a numerical simulation of the 
RG flow of the system of Eqs. [3] to [5] One can see that 
by taking 7V 7 >> Nz >> Ng we can achieve rj£ << r\ v . 
Because the RG flow of Zi v z is driven only by having a 
small number of fcrmionic species, the coupling constant 
runs only logarithmically. However, the huge number 
of gauge fields that participate in the RG equation of 
speeds for the charged particles force the running of r\ to 
be extremely fast. 

Although in the analysis conducted above we consid- 
ered a constant number of hidden particles, this number 
may depend on the energy scale. We can model this de- 
pendence as a power law 
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FIG. 2: Numerical simulation of the running of the charge 
Ztv,z, and ratios gi and g v . We use the initial conditions 
c * = c e — 2, c 7 = cz = 1 and Z* tvZ = Z^ ln = 1 at /i* = 100. 
We take JV 7 = 5 x 10 s , N z = 2 x 10 5 , N t '= 50, and { gi } = 1. 
We see that the running of Zi v ,z is logarithmic, while the 
running of gt and g v is power-law. Very small values of g are 
achieved in a very short interval of running. Moreover, we see 
the hierarchical structure of the final speeds: g e << g v . 
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FIG. 3: Numerical simulation of the running of the charge 
Zi v ,z, and ratios gi and g v using the power law dependence 
in Eq. [15] We take the IR mass values M e = M v = M 7 = 1, 



and we use I\ 



0.02, T z 



0.6, r. 



.5, and a 7 = 2.1, az 



0.65, on = 0.1. We use the initial conditions c* v 



2. 



c* — c* z — 1 and Z\ vZ = Z* Url = 1 at fj,* = 10 , and we set 
{gi} = 1. We note that the value of ge is saturated by the 
error tolerance of the code. More powerful computations may 
give smaller values. 



where Ti and a, are positive constants, and M$ is an 
IR mass scale. This exact behavior is also exhibited in 
models of large extra-dimensions where the Kaluza-Klein 
(KK) modes (from the 4D point of view) obey a power 
law as in Eq. [TS] [ll[. 1 In this context, M is the lowest 
KK mode Mi ~ 1/L, where L is the size of the extra 
dimension, and a — d is the number of extra dimensions. 



It is also interesting to note the equivalence between theo- 
ries with extra dimensions and strongly coupled QFT via the 
AdS/CFT correspondence. 
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In FIG. [3] we show the simulation of the RG system using 
the power-law formula in Eq. [15] 

To understand the choice of parameters used in the 
simulation in FIG. [3] it is instructive to calculate the 
total number of fermions and gauge fields as seen in the 
UV. Using Eq. [15] and the numerical coefficients given in 
FIG. Owe find A 7 - 10 11 , N z ~ 10 4 and Ne ~ 1, which 
explains the hierarchical structure of the final speeds. 

Now we comment on the validity of the RG equa- 
tions in a perturbative treatment of QFT. It is inevitable 
that one must make sure higher loop corrections are sup- 
pressed, otherwise our perturbative treatment (i.e. the 
loop-expansion) breaks down. As was stressed in [l[ , the 
condition for the validity of perturbation theory puts a 
severe constraint on the number of species. This con- 
straint may be relaxed if we notice that the condition for 
validity of perturbation theory is important mostly at the 
initial running, when the charges i?« )7 and Z^z are rela- 
tively large. However, once the charges run down the low- 
energy scale, for example by means of non-perturbative 
RG treatment, the condition for validity of perturba- 
tion theory may be satisfied. Thus, as long as the non- 
perturbative result does not change the overall trend dra- 
matically, fast running would be expected to still be ob- 
tained even if the perturbative analysis is not reliable in 
detail. 



4. DISCUSSION 

If correct, the OPERA result would obviously have a 
dramatic effect on the Lorentz symmetry of the funda- 
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